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Abstract 

We studied the dynamical influence of a dust component on the gaseous phase in central regions of galactic 
disks. Therefore, we performed two-dimensional hydrodynamical simulations for flat multi-component disks 
embedded in a stellar and dark matter potential. The pressure-free dust component is coupled to the gas 
by a drag force depending on their velocity difference. 

It turned out that the most unstable regions are those with either a low or near to minimum Toomre 
parameter or with rigid rotation, i.e. the central area. In that regions the dust-free disks become most 
unstable for a small range of high azimuthal modes (m ~ 8), whereas in dusty disks all modes have similar 
amplitudes resulting in a patchy appearance. The structures in the dust have a larger contrast between 
arm and inter-arm regions than those of the gas. The dust peaks are frequently correlated with peaks of 
the gas distribution, but they do not necessarily coincide with them. This leads to a large scatter in the 
dust-to-gas ratios. The appearance of the dust is more cellular (i.e. sometimes connecting different spiral 
features), whereas the gas is organized in a multi-armed spiral structure. 

We found that an admixture of 2% dust (relative to the mass of the gas) destabilizes gaseous disks 
substantially, whereas dust-to-gas ratios below 1% have no influence on the evolution of the gaseous disk. 
For a high dust-to-gas ratio of 10% the instabilities reach the saturation level already after 30 Myr. 
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1 Introduction 

Chemodynamical simulations are characterized by a multi- 
phase treatment of the interstellar medium (ISM). In most 
of these calculations the ISM is split into two dynamically 
different components, the clumpy star-forming molecular 
clouds and the diffuse warm or hot gas (Theis et al. 1992, 
Samland et al. 1997). More refined decompositions of the 
ISM are only done for one-zone models (e.g. Ikeuchi et 
al. 1984), but not for gravitationally coupled, dynamically 
evolving systems. Especially, the dust component in galax- 
ies is usually just seen as a tracer population which has no 
dynamical influence on galactic dynamics. On the other 
hand, it is well known that even a small admixture of a 
dynamically cold component can substantially destabilize 
galactic structures as Jog & Solomon (1984) demonstrated 
for a small amount of cold gas embedded in stellar disks. 
The aim of our paper is to study, whether a dynamically 
cold dust phase can also destabilize galactic disks. 

Good places for looking on the influence of dust are 
galactic central regions. During the last years high resolu- 
tion observations revealed that the circumnuclear regions 
contain stellar-gaseous mini-disks with sizes of a few hun- 
dred pc (e.g. CaroUo et al. 1997). These disks have an as- 
tonishing rich structure described as mini-bars, spiral-like 
dust lanes, star-forming rings or spiral arms. Many nuclear 
regions are also well described as patchy or multi-armed. 
Investigations by Regan & Mulchaey (1999) showed that 
these spirals are the most common morphological struc- 
tures in the central regions of galaxies. 

Since the nature of the mini-spirals is of great impor- 
tance for our understanding of a variety of astrophysical 
processes like the mass accretion into galactic nuclei, sev- 
eral explanations for have been invoked. E.g. Athanassoula 



(1992) studied the gas flows in and around bars by (stellar- 
dynamical) orbital analysis. She found that the existence 
of different periodic orbit families (or equivalently the exis- 
tence of inner Lindblad resonances (ILR)) is a key criterion 
for the existence and the shape of dust lanes. 

Different to the stellar-dynamical interpretation Engl- 
maier & Shlosman (2000) suggested that mini-spirals in 
central regions of galactic disks are related to the forma- 
tion of grand-design spiral patterns in galaxies. They ar- 
gued that gas density waves - different to stellar density 
waves - are not completely damped or absorbed at the ILR 
and, thus, they may generate spiral structures at all radii 
including the nuclear regions. Such a model might explain 
the continuity of some spiral features at small and large 
radii as well as the low arm-interarm contrast observed in 
galactic centers. A similar idea of induced structure for- 
mation is to invoke secondary bars or small nuclear bars 
located inside the ILR (Wada & Koda 2001). 

A common property of all the mentioned mechanisms 
is that they result in structures dominated by two or a few 
arms. On the other hand, Elmegreen et al. (2002) stressed 
that nuclear dust spirals differ from main-disk spirals in 
several respects: the nuclear spirals are very irregular with 
both trailing and leading components that often cross. 

All discussed mechanisms have some weak points and 
can not fully account for all the observational data. We 
consider here a new approach based on the observation that 
circumnuclear disks of galaxies are dusty. By means of a 
stability analysis Noh et al. (1991) showed that a dust com- 
ponent can strongly destabilize proto-planetary disks. The 
admixture of only 2% of dust enhances the growth rates of 
the dominant gaseous phase signiflcantly. A conservative 
estimate of the dust-to-gas ratio r in the solar neighborhood 
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Figure 1: Perturbations of the surface density (normalized to their initial values) of the gas and dust component 
after t = 6 ~ 90 Myr for a model with a 10% dust mass fraction (relative to gas): ASg^^/Sg^d for gas (left) and 
dust (right). Areas devoid of material are white, areas with a density enhancement of a factor of 2 or more are 
black. The grey area seen at the outer edges correspond to no deviation from the initial surface density. The 
lengths at the x— and y— axes are given in kpc. 



gives an average value of 0.6% ranging from 0.2% up to 4% 
in H2 regions (Spitzer 1978). These are lower limits because 
large, massive dust grains do not have a detectable extinc- 
tion. Direct satellite-based measurements give r-values up 
to 2% (EYisch et al. 1999). Little is yet known about the 
true dust-to-gas ratio in external galaxies, but some galax- 
ies show large values. E.g. M51 has a dust-to-gas ratio of 
2%. Taking the observed relation between the dust-to-gas 
ratio and the metallicity (Issa ct al. 1990) into account, it 
is reasonable to assume larger dust mass fractions in the 
metal-enriched central regions of galaxies than on average. 

In this paper we study the dynamical influence of a 
dust component on circumnuclear gaseous disks. We re- 
strict our analysis to systems without a nuclear or a large 
scale bar and we do not consider any formation or destruc- 
tion processes of the dust. This investigation is done by 2D 
multi-component hydrodynamical simulations. 



2 Method 

2.1 Equations of motion 

We studied numerically the hydrodynamical equations for 
a 2-dimensional single- or multi-component disk. Thus, we 
solved the continuity equation 

^ + V-(E,,di7,,d)=0 (1) 

and the momentum equations for gas and dust. The mo- 
mentum equations read for gas 

riv VP 

^ + {^9- V)iT« + ^ + V($ + *hbsd) = S,{v,) (2) 
and for dust 

^ + {^d- V)iJd + V(* + 3.HBSD) = Sd{Vi) . (3) 

The components are denoted by g for the gaseous phase 
and by d for the dust. Eg^^ are the surface densities and 



Vg^d the velocities. Pg is the pressure of the gas which is 
given in the 2d-case as force per unit length. $ denotes the 
potential of the self-gravitating disk. It is derived from the 
Poisson equation 

A* = 47rGE(i?, ifi) S{z) = 4nG (E^ + Ed) 6{z) . (4) 

An external potential is added by a stationary contribution 
$HBSD related to the halo, the bulge and/or a stellar disk 
component. These external potentials are chosen to match 
- together with the potential of the disk - a given rotation 
curve. 

The main difference between the gaseous and the dust 
component is the treatment of the dust as a prcssureless 
phase. Therefore, if gas and dust are in rotational equi- 
librium, there is a velocity difiierence between both com- 
ponents which might give rise to a non-negligible frictional 
force depending on the cross-section for the gas-dust inter- 
action. This interaction is described by the source terms 
S{. . .) on the RHS of the hydrodynamical equations. Since 
we do not consider dust formation and destruction pro- 
cesses, the source terms in the continuity equation vanish. 
However, frictional terms show up in the equations of mo- 
tion. The dust implementation will be described in the next 
paragraph. 

The set of hydrodynamical equations is closed by a poly- 
tropic equation of state 

Pg^KT.-'g^ . (5) 

For the gaseous phase we set the polytropic exponent to 
fg = 5/3. The constant K is chosen to yield a given mini- 
mum Toomre parameter. 

2.2 Treatment of the dust component 

The main difference between "normal" galactic disks and 
our dusty disks is that the cold component (dust) is not 
only coupled by gravity to the hotter phase (gas), but also 
by a frictional force between both components. This drag 
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is taken into account by a source term in the equations of 
motion following the general form suggested by Noh et al. 
(1991) 



'A{vd - Vg) 



(6) 
(7) 



The second source term, Eq. (7), follows from the require- 
ment of momentum conservation. The physics of the fric- 
tion is enclosed in the frictional timescale . Based on a 
microscopic view of the frictional process, i.e. the momen- 
tum exchange between gas and dust particles by collisions 
and the equipartition of momentum within the gaseous 
phase, the timescales can be calculated. For the simulations 
shown in this paper we assumed that the dust disk is thin 
compared to the gaseous disk. In that case the frictional 
time scale is of the order of the dynamical time scale given 
by the inverse circular frequency Sl^^ (Noh et al. 1991), i.e. 



A = Q. 



(8) 



The basic assumption is here that gas and dust establish 
very fast coUisional equilibrium in the thin dust layer and 
that the gas momentum is then mixed vertically in a sound 
travelling time scale (for further details see Theis & Orlova 
2004). 



2.3 Numerical implementation 

The nonlinear analysis implies the solution of the full set of 
hydrodynamical equations. For this purpose, we developed 
a two-dimensional numerical code which is similar to the 
ZEUS-2D code by Stone & Norman (1992). The hydrody- 
namical equations are discretized on a logarithmic Eulerian 
grid in polar coordinates (270x270 grid cells). The differ- 
ent terms are treated by operator splitting. Advection is 
performed by a second order Van Leer advection scheme. 



3 Results 

3.1 Parameters of the initial gaseous disk 

The characteristic values of our initial models are moti- 
vated by the nuclear region of MlOO (NGC 4321). We 
adopted a total gas mass of 4.7 • 10* M© distributed expo- 
nentially with a scale length of 300 pc within a radial range 
of -Rin = 30pc to Tiout = 3kpc. The short disk scale length 
mimicks a central concentration of (molecular) gas observed 
in many galaxies. Small random perturbations are initially 
superimposed on the mass distribution by multiplying the 
surface density in each cell with a factor (1 + ArR). R 
is a random number in the interval [-1,1] and Ar a small 
amplitude of the order of 10~*. 

The initial velocities are derived from rotational equi- 
librium, i.e. there is no initial radial motion. The rotation 
curve is specified by 



Vc{R) = «o 
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in a fairly sharp transition at the radius -Rflat = 100 pc. 
This rotation curve corresponds to a total dynamical mass 
Md(R) ~ vl{R)R/G (including aU components) of 7.3 • IO^Mq 
within the central 100 pc. In the region of rigid rotation 
the rotation period is about 3.5 ■ 10® yrs. It increases out- 
wards reaching 1.7 • 10^ yrs at the half-mass radius of the 
gaseous component at ii ~ 500 pc. The azimuthal veloc- 
ity of the gaseous phase is calculated by the (frictionless) 
Jeans' equation. 

According to the chosen rotation curve, mass profile 
and equation of state (polytropic with 7^ = 5/3), the mini- 
mum value of the Toomre parameter is reached at a galac- 
tocentric distance of about 440 pc. The constant in the 
equation of state (for the gaseous phase) was selected to 
yield a minimum Toomre parameter of Qmin = 1.54. This 
choice corresponds to sound speeds between 4 and 11 km 
s~^ within the central kpc (the higher value is reached in 
the center). 
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The velocity at infinity, Vaa, was set to 178 km s . The 
transition parameter rit was selected to be 10 resulting 



Figure 2: Temporal evolution of the logarithmic 
Fourier amplitudes Cg = log of the (dominant) 
m = 8-mode of the dust component for different dust- 
to-gas mass fractions: 0.5% (open boxes), 1% (crosses), 
2% (reference model, sohd Hue), 10% (plus), 20% (tri- 
angle) and the purely gaseous model (dashed line). The 
time unit is 1.5 • lO*" yrs. 



3.2 A dusty disk 

We studied the influence of the dust-to-gas mass ratio r 
by varying r = Md/Mg from 0.5% to 20%. Fig. 1 shows 
the surface density perturbations of gas and dust in the 
saturation stage of a model with a large dust-to-gas ratio: 
the irregular multi-armed structure is emphasized due to 
large arm-interarm variations. Along the arms strong sur- 
face density variations can be discerned. Some arms are 
interrupted by low density areas, others are not smoothly 
curved, but show wiggles. Some arms seem to merge with 
others. The dust distribution is highly correlated with the 
gas distribution. However, the contrast between arm and 
interarm regions is larger for the dust than for the gas. 
The surface densities of the dust vary by about one order 
of magnitude, whereas the contrast of the gas component is 
usually less than a factor of 2. The structures formed in the 
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dust are also thinner than those of the gas. Though there 
is a tight correlation between the positions of the maxima 
of the gas and dust phases, there is no correlation between 
their maximum amplitudes. The dust distribution is char- 
acterized by a more cellular appearance compared to the 
spiral-like morphology of the gas. 

It is interesting to note that the dust is often located at 
the boundaries around peaks of the gas distribution, pref- 
erentially at the inner boundary. Other dust peaks are in 
regions with no or only weak gaseous density enhancements. 
And, of course, there are dust peaks at the same locations 
as those of the gas. Thus, one expects large spatial varia- 
tions of the dust-to-gas ratio. 

In order to quantify the (global) stability of the disk we 
use the global Fourier amplitudes of each component 

Am = j^ I I E(r,0)rdre-"""^ci0 (m > 0). 



(10) 

Afdisk is the mass of the disk for the component of interest 
in the specified radial interval [flin,-Rout]- S denotes the 
corresponding surface density. Comparing the modes from 
m = 1 to m = 16, the high-m-modes around m ~ 8 are 
dominant during the linear regime. A comparison of the 
corresponding Fourier amplitudes with those of a purely 
gaseous model shows that the critical dust-to-gas ratio Vc 
for dust becoming dynamically unimportant is about 1% 
(Fig. 2). For larger amounts of dust the dcstabilization of 
the gaseous phase becomes much stronger. For r = 0.02 the 
instability sets in about 50 Myrs earlier than in dustless or 
low-dust-content models while the modes remain on their 
initial value for a latency period of about 75 Myrs. Increas- 
ing r to 10% reduces the latency time to almost zero. The 
growth rates increase by a factor of 3-4 and the saturation 
level reached already after 30 Myrs is larger by at least one 
order of magnitude. 

A more detailed discussion on the method and the re- 
sults as well as an extended parameter study can be found 
in our forthcoming paper (Theis & Orlova 2004). 



Summary 

We studied the influence of a cold dust component on the 
evolution of galactic gaseous disks by means of 2-dimensional 
hydrodynamical simulations. Prom the evolution of the 
Fourier amplitudes we found that the higher-order modes 

arc the dominant unstable modes. Their growth is mainly 
restricted to the central kpc. An admixture of 2% dust (rel- 
ative to the gas mass) destabilizes the gaseous disk. The 
growth of instabilities in the dust component becomes non- 
linear after 100 Myr, followed by the gas 50 Myrs later. 

The structures formed in both, gas and dust, arc rather 
irregular and multi-armed. This is a direct consequence 
of the superposition of many high-m modes of similar am- 
plitude. The dust component is characterized by a much 
stronger contrast between arm and inter-arm regions than 
the gas. The dust is spatially correlated with gas, but it 
does not exactly follow the gas. Peaks in the dust distri- 
bution are frequently found at the inner edges of peaks in 
the gas distribution. This results in a large scatter of dust- 
to-gas ratios at different places. The dust develops also 



thin filaments which sometimes connect the arms. There- 
fore, the dust distribution has a more cellular appearance, 
whereas the gas develops a multi-armed spiral morphology. 

Below a dust-to-gas mass ratio of 1% the dynamical 
influence of the dust on the gaseous disk becomes negligi- 
ble. This critical value is close to the observed mean value 
in normal galaxies like the Milky Way. Since the dust-to- 
gas ratio scales linearily with metallicity, larger local values 
of r, especially in the central galactic regions, seem to be 
reasonable. Such values arc also in agreement with local 
gas-to-dust determinations. Since already a dust-to-gas ra- 
tio of 2% significantly affects the evolution of the disk, even 
the observed small dust admixtures are expected to have an 
impact on the dynamics of some galaxies (e.g. the dust-rich 
M51). For a 10% admixture of dust the gaseous component 
is completely destabilized. The growth rates are enhanced 
by a factor of 3-4 with almost no latency phase. The sat- 
uration levels reached after 30 Myr are substantially larger 
than in the low-r calculations. 
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